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Abstract 

The purpose of this paper is to study the relations between different 
concepts of dispersive solution for the Vlasov-Poisson system in the grav- 
itational case. Moreover we give necessary conditions for the existence 
of partially and totally dispersive solutions and a sufficient condition for 
the occurence of statistical dispersion. These conditions take the form of 
inequalities involving the energy, the mass and the momentum of the solu- 
tion. Examples of dispersive and non-dispersive solutions — steady states, 
periodic solutions and virialized solutions — are also considered. 

AMS classification (2000). Primary: 35B05, 35B40, 82B40, 82C40 

1 Introduction 

Gravitational systems composed by a large number of classical particles in which 
collisions and external forces are negligible, e.g. the stars of a galaxy, can be 
described by the Vlasov-Poisson system. The regularity of solutions to this 
system has been extensively studied in the mathematical literature and this 
problem is by now well-understood, cf. [IHl HH 131] ■ On the contrary, very little 
is known on the time asymptotics of the Vlasov-Poisson system. 

The large time behavior of solutions to the Vlasov-Poisson system is rela- 
tively simple in the case of electrostatic interaction among the particles, which is 
obtained from the gravitational case by reversing the sign in the right hand side 
of the field (Poisson) equation. In the electrostastic case all solutions exhibit 
a (strong or L'^-norm) dispersive character [15l [21] . In the gravitational case 
the dynamics is more intricate: There exist (stable or unstable) steady states, 
periodic solutions (breathers) and (fully or partially) dispersive solutions. For 
the applications in Astrophysics it would be desirable to have a classification 
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of the possible asymptotic behavior of solutions in terms of relations between 
quantities preserved by the evolution (such us the energy and the mass). This 
is clearly a very difficult — may be impossible — task, but in this paper we show 
that partial answers in this direction can be given for the Vlasov-Poisson system. 
We shall focus in particular on solutions that for large times exhibit some sort 
of dispersive behavior. The examples considered in the following sections show 
that solutions of the Vlasov-Poisson system may disperse in several different 
ways. Another goal of the present paper is to study the relations between these 
various concepts of dispersion. 

Let us recall some fundamental facts about the Vlasov-Poisson system in the 
gravitational (galactic dynamics) case (more details can be found in [H]). The 
dynamics of the stars of the galaxy is described by the distribution function in 
phase space / = f{t,x,p), which gives the probability density to find a particle 
(star) at time t G R in the position a; e M"^ with momentum p £ . The mass 
density p = p(t, x) of the galaxy is given by 

fdp. (1.1) 

For notational simplicity we assume that the stars have all the same mass m 
and fixed units such that m — AttG = 1, where G is Newton's gravitational 
constant. The gravitational potential U = U{t,x) generated by the galaxy 
solves the Poisson equation 

A^U^p, lim [/ = 0,VteR, (1.2) 

\x\ — ^oo 

where the boundary condition at infinity means that the galaxy is isolated. The 
assumption that the stars interact only by gravity leads to the Vlasov equation: 

dtf+p-y.f -y.u -Vpf ^0. (1.3) 

The system (|l.ip - (|1.3p is the Vlasov-Poisson system. The solution of (|1.2p is 
given by the formula 

U = -^f ^,dy, (1.4) 
47r \x - y\ 

whence the Vlasov-Poisson system is equivalent to the non-linear Vlasov equa- 
tion obtained by replacing the formula for U in (|1.3p . The energy E and the 
mass M of a solution are given by 

E=l [ \p\\fdpdx-l [ \V^U\^dx , M^f fdpdx (1.5) 

2 JR6 2 Jjj3 

and are conserved quantities. Likewise, the total linear momentum Q and an- 
gular momentum L, 

Q — p f dpdx , E — X Ap f dpdx , (1-6) 
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are conserved quantities. Moreover, since the characteristic flow of the Vlasov 
equation preserves the Lebesgue measure, then all L'' norms of / are preserved: 

= const. , for all 1 < q < cx) . (1.7) 

Hereafter we denote by either L^(R^) or L^{M?), depending on the function 
under consideration, e.g. f{t) = fit., ■, ■) or p{t) = p{t, •). The invariance of 
Vlasov-Poisson by (time dependent) Galilean transformations is the property 
that, given u G M"^ and the transformation of coordinates 

Qu '. t' ^ t , x' = X — ut , p = p — u , (1-8) 

then fu{t,x,p) = f{t',x',p') and Uu{t,x) — U{t',x') solve the system (|l.ip - 
(|1.3p if and only if (/, U) does. Note that Q can be made to vanish with a 
suitable Galilean transformation; the resulting reference frame is at rest with 
respect to the center of mass of the distribution, which is defined as 

Cp{t) = M-^ / xpdx . (1.9) 

Throughout the paper we assume that / is a non-trivial global classical solution 
of the Vlasov-Poisson system such that, at any fixed time t, f has compact 
support in the variables {x,p) (however this assumption can be substituted by 
suitable decay conditions on the variables (x,p) or by requiring only that / has 
bounded moments in these variables up to a sufficientely high order). We shall 
refer to these solutions as regular solutions. It is well known that for any initial 
datum < — f\t=o G C^ci^^)-, there exists a unique global regular solution 
of the Vlasov-Poisson system, sec [20, 22j ^31^ . 

An important role in our discussion is played by spherically symmetric solu- 
tions of the Vlasov-Poisson system. A solution of Vlasov-Poisson is spherically 
symmetric if f{t, Ax, Ap) — f{t, x,p) for all rotations A G S0(3). The potential 
induced by a spherically symmetric solution is a function of the radial variable 
r = |a;| only and indeed we have the representation formula 

1 r 

drU =— / X^p{t,X)dX. (1.10) 
^ Jo 

Clearly, the center of mass of spherically symmetric solutions is at r = 0. 

This paper is organized as follows. In Section[2]we introduce several concepts 
of dispersion for a mass distribution — not necessarily originated by the Vlasov- 
Poisson system and give some examples. In Section [3] we specialize to the case 
of the Vlasov-Poisson system. We give necessary or sufficient conditions for 
the existence of various types of dispersive solution. For completeness we also 
briefly discuss in Section U] two examples of non-dispersive solutions, namely 
periodic solutions and steady states. 

To conclude this Introduction we remark that the Vlasov-Poisson system 
ceases to be valid as a physical model when the stars move with large velocities 
(of the order of the speed of light) or in the presence of very massive galax- 
ies, since then relativistic effects become important. Typical relativistic effects 
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are the redshift of the luminous signals emitted by the galaxy and the forma- 
tion of black holes. The model which is currentely believed to represent the 
physically correct relativistic generalization of the Vlasov-Poisson system is the 
Einstein- Vlasov system [I] , where Poisson's equation is substituted by Einstein's 
equations of General Relativity. As compared to the Vlasov-Poisson system, the 
Einstein- Vlasov system is far more complicated and less understood. Another 
relativistic generalization of Vlasov-Poisson is the Nordstrom- Vlasov system [7] , 
which, although not physically correct, is mathematically interesting since it al- 
ready captures some of the technical and conceptual new difficulties that are 
encountered when studying a relativistic {i.e., Lorentz invariant) system. The 
large time behavior of the relativistic models will be discussed elsewhere [5] . 

2 Dispersive behavior 

In this section we introduce several concepts of dispersion for regular mass 
distributions. By a regular mass distribution of total mass M we mean a non- 
negative function p{t, x) such that p{t, •) has compact support and ||p(i)||i = 
M (independent of time t). This terminology is consistent with the one used for 
solutions of the Vlasov-Poisson system: The mass distribution p defined by 
is regular whenever / is regular. 

2.1 Strong dispersion 

Definition 1. A regular mass distribution p is said to be strongly dispersive if 
there exists q > 1 such that the limit 

lini ||p(t)||L<! exists and is zero. (2-11) 

t — >oo 

Obviously, strong dispersion is a Galilean invariant concept. For the Vlasov- 
Poisson system in the plasma physics case, which is obtained from (|l.ip - (|1.3p 
by reversing the sign in the right hand side of (|1.2p . it was proved in [121 [H] (see 
also [TT]) that all solutions are strongly dispersive, with (|2.1ip being verified for 
q G (1,5/3]. In the gravitational case, examples of strongly dispersive solutions 
are those constructed in [3] for small initial data, see also [Mj . For these solutions 
there holds the estimate 

p{t,x)<C{l + \t\ + \x\)-\ 
for a positive constant C, which clearly implies strong dispersion. 

2.2 Total and partial dispersion 

The next types of dispersive solution that we are going to discuss use the notion 
of 'concentration function of a measure' introduced by Levy ^17j and applied by 
P.-L. Lions in the proof of the concentration-compactness Lemma |19j . 
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Definition 2. A regular mass distribution p is said to he totally, respectively 
partially dispersive, if and only the limit 

M{R) = lim sup / p{t,x)dx, (2.12) 

exists and 

Moo = lim M{R) (2.13) 

R — ^oo 

satisfies Aioo — 0, respectively A4oo £ (0, M). 

Remark 1. Of course, it is possible that Ai{R) could not be well defined for 
all R {e.g. when p{t) is time periodic). Whenever it exists, A4{R) is a bounded 
non-decreasing function and therefore the limit (|2.13|) is well defined. Moreover 
Moo G [0,M]. 

It is clear that strong dispersion implies total dispersion. Moreover, total dis- 
persion is equivalent to the vanishing property in the concentration-compactness 
theory, see [TH]; precisely, a mass distribution p is totally dispersive if and only 
if the limit 

lim sup / p{t, x) dx exists and is zero, Vi? > . (2-14) 

*^°°xoeR=' J\x-xo\<R 

An important physical property of (I2.14p is that it is invariant by Galilean trans- 
formations, unlike the decay of the mass (or energy) over a ball with arbitrary 
radius, 

/ pdx^O, ast-^oo, Vi?>0, (2.15) 

J\x\<R 

which has also been used as definition of dispersion for evolution type equations 
(including non-linear Vlasov equations), see [11|, I33j for instance. Thus for 
example, according to our definition of total dispersion, a static (i.e., time 
independent) solution which is "put in motion" by a Galilean transformation is 
not to be regarded as a dispersive solution (whereas it would be so according 
to (EUl)). 

Example 1. As an example of totally dispersive solution of the Vlasov-Poisson 
system, consider a spherically symmetric shell of matter with internal radius 
Ri(t) and — possibly infinite — external radius R2{t) (this example was first in- 
troduced in [2]). Let r = \x\ and w = x -p/r, the radial momentum variable. 
Now suppose that in the support of = f{0,x,p), it is verified that 



inf{.,.e Buppn>^^, (2.16) 

where M is the total mass, i.e., initially the particles are moving outwardly 
with sufficiently high speed. Using that in spherical symmetry the maximal 
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force experienced by a particle is bounded by M/47rr^, see (|1.10p . we find that, 
along the characteristics of the Vlasov equation, 

d / 1 2 Af \ M f M \ 

t: ~ — ] =ww + - — ^ r = w w + - — - , 

which is positive in the time interval [0,T) in which w > 0, i.e., as long as the 
shell keeps moving outwardly. It follows that 

w{tf > w{Of ^ > inf -W>0, 

^ ^ ^ ' 27rr(0) supp f 27ri?i(0) 

where for the last inequality we use (I2.16p . This implies that T ~ oo, that is, 
the shell moves outwardly for all future times. Moreover, > is a uniform 
lower bound on the radial momentum, which entails 

R2it)> Ri{t)> Ri{0) + Wt. (2.17) 

We claim that, because of (|2.17[) . the solution under consideration is totally 
dispersive. We shall achieve this by proving that the potential energy vanishes 
in the limit t oo, which for a solution of Vlasov-Poisson is equivalent to total 
dispersion, see Proposition [2] in the next section. Thanks to the rotational sym- 
metry we have the representation (|1.10p , which allows to estimate the potential 
energy as 

^^pot(i) = ^ / \V,U\^dx = i / f \^p{t,\) d\\ dx 

R2{t) ^ / r \ 2 



2tt ^ / \'^p{t, A) dX dr 



'fli(t) 

< — / — dr < 



and the claim follows. 

For the next result we denote by d{A, B) the distance of the sets A, i? C K^: 

d{A,B) = inf{|x - y\ ,x ^ A ,y ^ B} . 

Lemma 1. Let p he a partially dispersive regular mass distribution. Then, for 
any given e > there exist i„ — > oo and two sequences of non-negative 
functions p", p2 '-^^ ~* I^+i such that p{tn) > p" + p2 and 

b) |||p?||i-A^oo| <e,||lp^||i-(M-A^^)| <e; 

c) d(suppp", suppp2 ) ^ oo, as n ^ oo. 
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d) There exists a sequence of vectors G and < R*^^ such that — 0, 
for > R^y 

Remark 2. The conditions a)-c) define the dichotomy property of the mass 
distribution p in the concentration-compactness Lemma, see [19 . Condition 
d) is also consequence of the same result, although this was not pointed out 
in [19] ; for the sake of completeness we shall give here the proof of Lemma [TJ 
In our context, the relevance of the extra condition d) arises from the fact that 
it prevents the system from being strongly dispersive, as we will see in Section 

[sm 

Proof. The following proof is adapted from [T9^. Owing to (|2.14[) . along any 
sequence tn oo we have 

lim sup / p{tn, x) dx ~ M{R) . 

^^°°xom^ J\x-xo\<R 

Since M{R) Moo e (0,M), for aU e > we can find R* = R'^^^ such that, 
for all n sufficiently large. 



sup / p{tn,x)dx e {Moo - £,Moc + s) ■ 

KoSR^ J\x-xo\<R* 

Moreover there exists y" G M'^ such that 

p{tn, x) dx e {Moo - £, Moo + e) • 

x~y"\<R' 

Finally, we can find a sequence i?„ oo and a subsequence of times — still 
denoted tn — such that 

sup / p{tn,x)dx e {Moo - £,Moo + e) ; 

xo€R^ Jlx-xolKRn 

The functions p"^ = p{tn)x{By,^^R*)} and p'^ = p(tn)X{R3\s„,. (j?„)} are easily 
seen to satisfy the properties a)-d). □ 

2.3 Statistical dispersion 

We shall now discuss another Galilean invariant concept of dispersion, which 
was introduced in [10 . Define the statistical dispersion operator in space by 



/ \x\^p{t,x) dx - ^ ( ; 

JR3 



2' 

xp{t, x) dx 



Up to a mass normalization, the statistical dispersion operator coincides with 
the statistical variance of the density mass function and, consequently, it is a 
measure of the dispersion of such distribution. Note also that 

((Ax)^)= / \x-c,{tf {-dx 
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and therefore the statistical dispersion operator coincides with the moment of 
inertia of the mass distribution with respect to the center of mass. 

Definition 3. A regular mass distribution p is said to be statistically dispersive 
if and only if 

sup({Aa;)^) = +00 . 
t>o 

Remark 3. The above definition differs slightly from that given in [10], where 
statistical dispersion is defined by the condition limt^oo ((Ax)^) = +00. 

Statistical dispersion is the weakest concept of dispersion among those in- 
troduced so far. 

Proposition 1. If a regular mass distribution is totally or partially dispersive, 
then it is statistically dispersive. In particular, total dispersion implies that 
limt^oo((Aa;)2) = +00. 

Proof. We prove first that total dispersion implies limt^tx3((Ax)^) — +00. Fix 
R > arbitrarily and write 



M — pdx + pdx 

J\x-Cp(t)\<R J\x-Cp{t)\>R 

< sup / pdx + pdx . 

xoem.^ J\x-xo\<R. J\x-Cp(t}\>R 

Assume the solution is totally dispersive. Then by (|2.14p there exists to = to{R) 
such that, for all t > to, 

f ^ M 

sup / pax < — . 

xoeM.^ J\x-xo\<R 2 

Thus for t > to, 

((Axf)>|-/ pdx>^, 

J\x-Cp(t)\>R ^ 

which yields the claim. To prove that partial dispersion implies statistical dis- 
persion we use the dichotomy property of partially dispersive solutions, see 
Lemma [H Let i„,p",p2 satisfy the properties of Lemma [T] Then 



{{Axy){tn)> \X-Cp{tn)\^p'ldx+ \x - Cp{tn)\' P2 dx 

J SUppp^^ J SUppp2 

>d(Cp(t„),SUpppr)2||p^||l-Hd(c,(i„),SUppp?)2||p5||l 

> C [rf(cp(i„),supp/9'j')2 + d(cp(i„),suppp^)2] . 
By the triangle inequality, 

[. . . ] > (i(suppp", suppp27^)^ — > 00 , as f — > 00 , 
whence ((Ax)^)(t„) ^ oo, which concludes the proof. □ 
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Example 2. We give now an example of solution which is partially (and there- 
fore statistically) dispersive but not strongly dispersive. This example is a 
modification of the fully dispersive shell considered before, in which a static, 
spherically symmetric configuration with given mass Mq is located in the in- 
terior of a shell with mass m (alternatively, the interior part may consist of a 
static shell |25j . leaving a neighborhood of the origin empty, or a spherically 
symmetric periodic solution, such as the one found by Kurth |16| . see also next 
section). Since the potential inside the shell is constant, the static configuration 
in the interior will persist as long as the shell is moving outwardly. This again 
will be verified under condition (j2.16p . which now reads 



mf{w,w e supp /shell} > 

Then we have 

(Mo + m)((Aa;)2) / \x\'^pdx+ / {xl"^ pdx > Ri(tfm . 

By (|2.17p , this gives a growth of the spatial variance of order . Partial disper- 
sion also follows immediately by (|2.17p . 



i Mg+m 
V 27ri?i(0) ■ 



3 Dispersion in the Vlasov-Poisson system 

Recall that strong dispersion implies total dispersion. Indeed we shall now 
prove that these two concepts of dispersion are equivalent for the Vlasov-Poisson 
system. Let 

E^in = 11 \p\''fdpdx , Spot - ^ / \V^U\^dx 

and note that 

Skin-:^-^/ \p~M-'Q\'' fdpdx>0. (3.18) 

Proposition 2. Let f be a regular solution of the Vlasov-Poisson system. Then 
the following assertions are equivalent: 

1. f is strongly dispersive. 

2. f is totally dispersive. 

3. Epot —^0, as t ~^ oo. 

Moreover, if any of the above holds, then f satisfies the inequality 

E> — . 
- 2M 
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Proof. By (15351) . 



pot 5 



and the last claim follows by letting t ^ oo. Let us prove the equivalence of the 
three statements: 

1. =^ 2. Clear. 

2. 3. Fix R > and rewrite the potential energy as 87ri5pot = I1+I2+I3, 
where 

J J\x^y\<l/R \x-y\ 

— , , — dxdy, 

J Jl/R<\x~y\<R \X - y\ 

P^ifP^d.dy. 

J J\x-y\>R \x-y\ 

Using the Young inequality, see [TS], the first integral is bounded as 

we recall that ||/3(i)ll5/3 ^ for regular solutions of the Vlasov-Poisson system, 
see for instance [2l]- For J3 we use that 

^ - R 

Finally 

h < R p{t,x)p{t,y)dxdy < MR swp / p{t, x) dx — ReR{t) , 

J\x-y\<R ym'' J\x~y\<R^ 

where, by (|2.14p . £fl(t) — > 0, as t ^ 00, for all R> 0. Collecting, 



Taking the limit t ^ 00 and then i? — s- cx) concludes the proof. 

3. =^ 1. We recall from [15] the following interpolation inequality 



IIPII5/3 < Ct-^ (^^ \x-tp\'f(t,x,p)dxdi^ 



(3.19) 



and the pseudoconformal law for the attractive case 

4/ \x-tp\'^f{t,x,p)dxdp^ 
dt , „6 at 



^ \x-tp\^J{t,x,p)dxdp^ ^ ie\\V^U{t)\\l) ~t\\V,U{t)\\l. (3.20) 
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Integrating (|3.20p . we get 

/ \x-tp\^fdxdp- f Ixl""/" dx dp = t^\\W^U{t)\\l- f s\\\/,U{s)\\lds, 



so that 



0<t-^ [ \x-tp\^fdxdp<t-^ [ \x\^fdxdp+\\V^U{t)\ 



2 

2 J 



and the r.h.s. converges to zero by hypothesis, which in combination with (j3.19p 
concludes the proof. □ 

In [T^ a sufficient condition for statistical dispersion was estabhshed, which 
we reprove here in a slightly different way. We shall use the identity 



M^(iAx)^)=2E + 2E^,^-2^ 
dt'^ M 



(3.21) 



which is proved by direct calculation. 

Proposition 3. Solutions of Vlasov-Poisson which satisfy the condition E > 

2M 



O 

are statistically dispersive. In particular, there exists constants Ci,C2 > 
such that, for all sufficientely large times, 

Cit^ < ((Ax)2) < Cai^ 

Proof. First we rewrite p.2ip as 



(3.22) 



M^((Ax)2)=4£;-2|^ + 2i?pot, 
Using that -Epot > and integrating in time twice we get 

{{Axfm > {{Axfm 

where 



(3.23) 



t=0 



dt 



({Axf 



2 

M 



x-{p-M-^Q)f^{x,p) dxdp. 



(3.24) 



The bound from below follows immediately. To prove the upper bound, we 
recall — see [TOl HI] for instance — that the potential energy satisfies the bound 



in 1 

where the positive constant C depends only on M — \\f{t)\\i and ||/(^)||c 
||/°||oo. Thus 

E\dn — C\/ i?kin ~ E <^ , 
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which in the case of non-negative total energy E gives immediately a uniform 
upper bound on the kinetic energy: 



1 \ ^ 

-{c+Vc^TaeU (3.25) 



and therefore the potential energy is uniformly bounded as well. Using this 
in p.2ip it follows that ((Ax)^) < 6*2^^ and the proof is complete. □ 

The threshold / (2M) represents the kinetic energy of a point at the center 
of mass having the same mass of the whole system. The case E = Q'^/{2M) is 
settled in the following proposition. 

Proposition 4. Solutions of Vlasov-Poisson which satisfy E ~ are statis- 
tically dispersive. 

Proof. After a Galilean transformation we may assume (5 = 0. In this frame, 
the solutions under consideration have zero energy. Thus we compute 

^-^((Ax)2)= / {x-p)fUxdp+ f Eun{T)dT. 
^ "I JR6 Jo 

The claim is obvious if /(a; • p)f'^ dxdp > 0. More exactly, in the latter case 
we have ((Ax)^) > Ct, for a positive constant C. Otherwise there exist two 
possibilities: 

1. There exists t* > such that — J{x ■p)f^dxdp ~ /J E]^ia{T)dT. In this 
case too {(Ax)^) >Ct. 

2. There holds the bound Jg°° £'kin(T) c?t < — J{x ■ v)f^ dxdv. Hence there 
exists a sequence of times t„ oo such that -Ekin(^?i) 0, and thus 
Epot{tn) ^ as well. From the proof of the implication 3 =^ 1 in Propo- 
sition[2l we see that ||/o(t„) II5/3 0. In particular 



lim sup / p{tn,x)dx^O, VR>0 



and repeating the argument at the beginning of the proof of Proposition[ll 
we conclude that ((Aa;)^)(i„) ^ +00. 

□ 

Example 3. At this point it is interesting to reconsider the example of the 
shell surrounding a static configuration introduced at the end of Section [2l Say 
that the inner part of the solution has mass Afo and that the surrounding shell 
has mass m and initial inner radius Ri. Our previous computations show that 
the escape threshold associated with this configuration is 



Mq + T? 

27ri?i 
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see (|2.16p . so that any particle with initial radial momentum greater than this 
threshold will escape towards infinity. We shall now prove that it is possible to 
obtain an escaping shell even when the total energy of the system is negative. 
Note that the total energy E consists of the energy Eq of the interior part 
plus the kinetic energy of the shell minus the potential energy of the shell (the 
interaction energy has negative sign and as a consequence the term contributing 
to the potential energy can be not consired). Neglecting the last negative term 
and estimating above the kinetic energy of the shell we get 

E < Eo + imsup , 

^ shell 

where the supremum is taken in the support of the shell at time t = 0. The 
interior part has energy Eq < 0, since it is static (cf. Proposition [5]); thus in 
order to have the whole shell escaping to infinity while the total energy of the 
system remains negative we must choose the initial radial momenta for all the 
particles in the shell according to 

Mo+m . r 2 • r I |2 I |2 -2i?0 



< inf < inf |pp < sup |pp < 

shni] «ho]I _i_„ii 



ZTTiti shell shell g^icW ^ 

This can be done if m is strictly contained in the interval between zero and the 
value i[— Mo + \/Mq — IQttEoRi]. For bigger values of m we have no guarantee 
that the total energy can be kept negative. 

The previous example shows us that: 

• There are solutions that are partially (therefore statistically) dispersive 
with E < Q'^/2M, so that there is no simple way to extend the results of 
Propositions [3] and [31 By Proposition [2l these solutions cannot be totally 
dispersive. 

• We know that spherically symmetric solutions with E > disperse sta- 
tistically with a dispersion rate of (equivalcntly their spatial support 
spreads with a velocity of order t) . We will see in Section 13.11 that for 
solutions with E — this needs not to be the case. On the other hand 
the previous example shows that there are solutions with E < which 
also statistically disperse with a rate t^. So there is no evident relation 
between the admisible rates of dispersion and the sign of the energy. 

• There is no lower limit for the fraction of total mass of the system that is 
lost to infinity for a partially dispersive system. 



Using these ideas we can also show that there exist dispersive solutions 
as "close" as desired to a stable steady state. More precisely, given a 
polytropic steady state (see [IS] for details and notation) with mass M, 
polytropic index /i and L^~^^^f^-iiOTm J, which is stable in the sense of 
[29| , we can find solutions as described above that are partially dispersive 
and that remain in the stability region of the polytrope. (This is not a 
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contradiction since the quantity of mass that is lost to infinity is almost 
negligible.) These solutions correspond to initial data of the following 
type. Starting form the given polytrope we construct — by scaling — a sec- 
ond polytrope having mass M — m and L^+^/^'-norm J—j, for to, j positive 
and small. Then we add an outer shell of mass to and L^+^/'^-norm less 
or equal to j. We choose m,j in order that all the computations in Ex- 
ample [3] remain valid and that the total energy of the solution is as close 
as desired to the energy of the original polytrope. Then we can invoke the 
stability criterium in |29j . 

3.1 Kurth's solution 

We conclude this section by considering an explicit class of spherically symmetric 
solutions to the Vlasov-Poisson system found by Kurth in Jj6j (see also [13], 
[24]). The main idea of this model is to find a solution whose associated density 
is of the form 

pit, x) = (47r/3)- V(i)"'x{|.|<0(t)} (3.26) 
where the function is interpreted as the radius of the system. This function 
(j) solves the following ODE 

030" -1-0 = 1, 

subject to the initial condition 0(0) = 1 . We get solutions to the Vlasov-Poisson 
system exhibiting different types of behavior depending on the prescribed value 
of 0'(O): 

• If 0'(O) = 0, the solution is static; 

• If < If/*' (0)1 < 1, the solution is periodic in time; 

• If |0'(O)| > 1, the solution is strongly dispersive. 

Let us relate the above classification in terms of 0' (0) with the values of the 
energy E. Note first that the associated distribution function / can be chosen 
to be spherically symmetric so that Q = 0. Doing so, the energy of Kurth's 
solutions is given by, see [T3] . 

i? =^((0')2 + 0-2- 20-1) = ^(0'(O)2-1). 

Moreover M = 1, which follows integrating (|3.26p . Thus 

• If _E = —3/5 (<^ 0'(O) = 0), the solution is static (this steady state is 
studied in [5]). 

• If -3/5 < i? < (^ < \(j)'{0)\ < 1), the solution is time periodic. 

• If i? > (<^ |0'(O)| > 1), the radius of the system goes to infinity and, 
by p.26p . p in L'^, for all g > 1, i.e., the solution is strongly (and 
therefore also totally) dispersive. When E = (<^ |</''(0)| = 1), we have 
{{Axf) - When > (4=» |0'(O)| > 1), we get {(Ax)^) - t^, in 
agreement with Proposition [3l 
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Let us prove the latter claim. First we show that ((Ax)^} ~ 0{(i){tY). Since 
the solution under study is spherically symmetric, we have 



((Ax)2} = / \x\^p{x) da; = 47r — 4>{t)-^ / dr ^ -(t>{tf ■ 

V 3 / Jq 5 

Thus we reduce the problem to find out the large time behavior of the function 
(j). Following Kurth J6J, if |(/)'(0)| = 1 we have that 



where v{t) solves 



v{t) + \v{tf ^2(1+"^ 



For t big the term v{tY dominates and then (t){t) ^ 0{t^/^). If > 1, we 

have that 

_ \<t>'mch{v{t))-i 

|<^'(0)|2-l ' 

where v(t) solves 

- |0'(O)|sh = (|</,'(0)| - - to) 

{to depends on |(/)'(0)|). For t big enough |sh(i;(t))| dominates and we 

infer that \v{t)\ ~ O(logi), which entails (j>{t) ^ 0{t). 

Note that the solution with E = (there are two of them actually) is the only 
known example of statistically dispersive solution for which statistical dispersion 
takes place at a slower rate than and for which the support spreads more 
slowly than t. Also the rate for strong dispersion is slower than for the other 
examples considered in this paper. Taking a closer look to the trajectories 
reveals that these are strongly oscillatory (like in a forced harmonic oscillator). 

This special solution highlights also the role of condition d) in Lemma [TJ 
For we can surround this "slowly dispersing" solution with a strongly dispersive 
shell configuration, and the resulting solution verifies a), b) and c) but not d), 
and happens to be totally but not partially dispersive, see previous Remark [2l 



4 Remarks on non-dispersive solutions 

In this section we briefiy consider non-dispersive solutions of the Vlasov-Poisson 
system. In particular, we consider steady state solutions and time periodic 
(breather) solutions. 

We distinguish between two types of steady states: Static solutions, which 
are defined as time-independent solutions of the Vlasov-Poisson systen (|l.ip - 
(|1.3|) . and travelling steady states, which are defined as solutions f{t, x,p) such 
that there exists a Galilean transformation Qu such that /ot/„ is a time indepen- 
dent solution of the Vlasov-Poisson system (i.e., a static solution). For static 
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solutions one has Q = 0, wheras Q ^ for travelling steady states. In fact, it 
is easy to show that the Galilean transformation which transforms a travelling 
steady state with linear momentum Q and mass M to a static solution (with 
linear momentum Q — and mass M) is the transformation Qu, where 



The construction of static solutions for the Vlasov-Poisson system can be 
done in two ways. Firstly, by choosing a particle density / which depends only 
on quantities that are conserved along the characteristics of the time indepen- 
dent Vlasov equation; with this choice the Vlasov equation is automatically 
satisfied and the problem is reduced to that of proving an existence theorem 
for the non-linear elliptic equation obtained by replacing / in the Poisson equa- 
tion. So far this 'direct' method was used mostly in the spherically symmetric 
case (see however [32]), where by the Jeans Theorem [5] all solutions of the 
time independent Vlasov equation can be expressed in terms of the particles 
energy and angular momentum. A second method to construct steady states 
of the Vlasov-Poisson system is by minimizing the energy (or a related) func- 
tional subject to suitable constraints (the choice of the functional and/or the 
constraints selects the type of steady state to be constructed). The advantage 
of the variational method on the direct method is that the former automatically 
proves a stability property for the steady state. We refer [T2J [551 [Ul UHl HI] and 
references therein for several works on the construction and stability of steady 
states of the Vlasov-Poisson system. 

To make a connection with the results proved in the previous sections for 
dispersive solutions, let us note the following simple bounds on the energy of 
steady states and time periodic solutions: 

Proposition 5. Static solutions of the Vlasov-Poisson system satisfy E < 0, 
whereas travelling steady states satisfy E < . Time periodic solutions of the 
Vlasov-Poisson system satisfy E < — 

Proof. The key to prove the result is the dilation identity: 



which follows by direct computation. If / is a static solution, then the previous 
identity implies the virial relation E = — i?kin, which yields the claim for static 
solutions. For travelling steady states, a Galilean transformation with u = 
(M)~-^(5 yields the virial identity in the form E — {M)~^Q^ = — i?kin, which, 
using (|3.18p . implies the claim of the proposition for travelling steady states. In 
the case of periodic solutions, we integrate the dilation identity over a period T 
to get 




(4.27) 




Using again ()3.18|) concludes the proof. 



□ 
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Remark 4. In the case of static solutions of the Nordstrom- Vlasov system, the 
inequality i? < is replaced hy E < M, see [8]. For the (spherically symmetric) 
Einstein- Vlasov system, one can obtain an inequality for static solutions which 
involves not only the energy and the mass but also the central redshift [5] ■ 

4.1 Virialized solutions 

It is a classical result in Astrophysics that bounded systems of self-gravitating 
particles roughly in equilibrium verify that the time average of the kinetic energy 
of the ensemble equals twice the time average of its potential energy. This 
statement and some of its variants and particularizations go under the name 
of "virial theorems" (cf. [23l [30l [6], for instance), and are common tools in 
Astrophysics. We shall comment here on the connection between the notion of 
virialized solutions of the Vlasov-Poisson system and our preceding results. 

In this section we shall only consider solutions such that Q ~ 0, i.e., the 
reference frame is chosen at rest with respect to the center of mass of the system, 
as we did in the proof of PropositionU) Following [23 we shall say that a solution 
of the Vlasov-Poisson system is virialized if and only if 

limii^±^^^^0. (4.28) 

t — ^oo t 

Note that solutions with E > Q cannot be virialized in this sense. In fact, it is 
a straightforward consequence of the inequality (I3.18P that virialized solutions 
of the Vlasov-Poisson system must necessarily satisfy 

E <Q . 

Examples of virialized solutions are Kurth's solutions with energy i? < 0. 

Remark 5. The notion of virialized system is usually applied in Astrophysics 
only in the case of iV-body bounded systems, but, as pointed out in [23], strict 
boundedness is not necessary (Kurth's solution with energy E = Q shows that 
also in the case of Vlasov-Poisson, the support of a virialized solution can spread 
out to infinity) . This remark led to interpret the virialized systems as apparently 
bounded systems, i.e., systems that disperse so slowly that in our time scale they 
appear to be bounded and in equilibrium. 

The following proposition extends the result in [23], which is valid for A^- 
body systems, to the continues setting; note that the diameter of the A^-body 
system used in [33], is replaced by the statistical dispersion operator ((Aa;)^). 

Lemma 2. Let f{t) be a given solution of the Vlasov-Poisson system with 
Q ~ 0. Then the following statements hold true: 

1. If f{t) is virialized then limf^oo ^^^^i ^ = 0. 

2. If limt^oo ^^^^i'^ ^ — and linit^oo -^o -E+-Ekin(T) dr ^^^^^^ then fit) is viri- 
alized. 



17 



Dispersive behavior 


Necessary 


Sufficient 


Example 


Strong Dispersion 

(VP)t i 

Total Dispersion 

i 

Statistical Dispersion 

T 

Partial Dispersion 


E > Q^/2M 
E > Q^/2M 


? 
? 


Example [T] 
Kurth's E>0 


7 


E > QV2M 




7 


? 


Example [2] 


Other solutions 








Static Solutions 
Periodic Solutions 
Virialized Solutions 


E <0 
E < -QV2M 
E <0 




Kurth's = -| 
Kurth's E 

Kurth's E = 



Table 1: Main results proved in the paper and open problems. 



Proof. The proof is a straightforward application of L'Hopital's rule (as formu- 
lated in Hg) and of fX^ . □ 

Remark 6. The existence of the limit in (|4.28p can be guaranteed under some 
of the most frequent situations: Spherical symmetry [4], periodicity in time 
(including static solutions of course) and whenever i?kin(i) has a limit as t — > oo. 

5 Summary and open problems 

The main results of this paper are summarized in Table [TJ Completing the 
entries marked with a question mark would lead to a considerable extension of 
the result presented here and of our understanding of the large time behavior 
of the Vlasov-Poisson system in the gravitational case. 
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